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ABSTRACT 


Two aspects of digital communication were investigated. In the first part, a 
FFT-based, M-ary FSK receiver in a Rician-fading channel was analyzed to determine 
the benefits of non-uniform windowing of sampled received data. When a frequency 
offset occurs, non-uniform windowing provided better FFT magnitude separation. 
The improved dynamic range was balanced against a loss in detectability due to 
signal attenuation. With large frequency offset, the improved magnitude separation 
outweighed the loss in detectability. An analysis was carried out to determine what 
frequency deviation is necessary for non-uniform windowing to out-perform uniform 
windowing in a slow Rician-fading channel. Having established typical values of prob- 
ability of bit errors, the second part of this thesis looked at improving throughput 
in a digital communications network by applying adaptive automatic repeat request 
(ARQ) protocols. The results of simulations of adaptive ARQ protocols with vari- 
able frame lengths is presented. By varying the frame length, improved throughput 


performance through all bit error rates was achieved. 
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I. INTRODUCTION 


Two aspects of digital communication were investigated. Chapters I, II] and IV 
present an analysis determining the advantages of non-uniform windowing in reducing 
the effects of frequency offsets caused by doppler shifts and other sources in a fading 
channel. Chapters V, VI and VII develop a simulation of adaptive automatic repeat 
request (ARQ) protocols to improve throughput performance through various bit 


error rates. 


A. EFFECTS OF NON-UNIFORM WINDOWING 

Digital communication involves the transmission of information in bit format. 
M-ary frequency-shift-keying (MFSK) uses M = 2* different frequencies to represent 
M different symbols (each containing k bits of information). A proper determination 
of the frequency sent enables the receiver to determine which symbol was transmitted. 

With the advent of real-time Fast Fourier Transform (FFT) processors, an 
MFSK receiver can be easily implemented. Selected output bins of the FFT will 
correspond to the M7 possible frequencies provided the sampling rate and the signal 
frequencies are closely related. In the case developed here, a bin separation of two 
will be used (proper determination of sampling rate and corresponding bin locations 
will be discussed in Chapter II). By comparing the magnitudes of the FFT output 
bins and choosing the bin with the largest magnitude, the transmitted symbol can be 
identified. 

Noise in the channel causes a variance in the magnitudes of the FFT output 
and an error can occur when the magnitude of one output bin is greater than the bin 


corresponding to the frequency of the signal sent. In the presence of frequency offsets 





(caused by doppler shifts or receiver oscillator drift, for example), the probability 
of error is further increased by a shifting of signal power into other frequency bins. 
Frequency offsets also contribute indirectly to erroneous frequency components by 
causing leakage. Leakage is encountered in FFT analysis and is a result of a periodic 
signal not being sampled over an integer number of cycles. A discontinuity will occur 
at the endpoints of the signal and will contribute erroneous frequency components to 
the FFT. These contributions can be seen in Figure 1.1(a) where the magnitudes of 


a FFT are given for a noiseless signal experiencing a frequency offset. 
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Figure 1.1: Effects of frequency offset with (a) rectangular, and (b) Hamming 
windowing. Dynamic range is improved using non-uniform windowing with a loss of 
signal detectability. 


Non-uniform windowing can be used on the sampled signal in order to smooth 
out discontinuities at the endpoints of the sampling interval and help reduce the effects 
of leakage. The reduced effects of leakage are balanced against signal attenuation 


caused by non-uniform windowing. With less energy, signal detectability deteriorates 


and the noise variance has a greater effect on proper signal decoding. This is seen in 





Figure 1.1(b) where a Hamming window is applied to the sampled signal. Although 


the dynamic range between the two components of interest has increased, signal 
strength has decreased. As the frequency offset becomes more dramatic, the benefits 
of reducing the amount of leakage by non-uniform windowing outweighs the loss in 
signal detectability. 

This thesis investigates the amount of frequency offset necessary to overcome 
the performance loss due to signal attentuation in a Rician-fading channel. The 
results of [1] are applied to a channel experiencing slow Rician fading to determine 
the effects of fading on the performance trade-off between non-uniform and uniform 
windowing. A statistical analysis is carried out to determine the probability of bit 
error as a function of signal-to-noise ratios as well as direct-to-fading ratios associated 
with Rician fading. From these results, the frequency offset at which non-uniform 
windowing out-performs uniform windowing is determined. It is shown that as the 
channel approaches Rayleigh fading (the direct-to-fading ratio decreases), the amount 
of frequency offset necessary to justify the use of non-uniform windowing decreases 


despite its signal attenuation characteristic. 


B. ADAPTIVE ARQ 


Having an understanding of the expected bit error rate experienced by the 
channel, Chapters V through VII investigate the possibility of improving throughput 
performance by implementing an adaptive automatic repeat request (ARQ) scheme 
in computer communication. 

Communication between computers can be described using the seven layer open 
systems interconnection (OSI) model given in Figure 1.2. Information from each node 
is first processed down to the network layer and into the data link control (DLC). 


The DLC then prepares the packet of information for transmission. Overhead bits 


are appended to the packet which include coding information for error detection and 
synchronization. In addition, a header and a tail is attached to the packet so that it 
may be identified by the proper receiving node. The processed packet, now a frame. 
is released from the DLC to the physical interface for transmission over the physical 
link. The transmitted frame is then received at the receiving node’s physical interface 


and a reverse process occurs. 


Data Link Control Data Link Control 


Physical Interface Physical Interface 


Physical Link 


Figure 1.2: In an OSI network communication representation, data packets are 
passed from the network layer to the data link control (DLC) layer where they are 
formatted and processed for transmission over the physical link. ARQ protocols are 
carried out in the DLC layer. 


As noted in the previous section, noise on the physical link sometimes corrupts 
the transmitted frame and an error-detection scheme is implemented by the DLC to 
recognize transmission errors. When an error occurs, the receiving DLC must issue a 


command to the sender of the frame to request retransmission of the frame. 





Several different ARQ protocols currently exist and three will be analyzed here— 
stop-and-wait, go-back-N and selective-repeat. An analysis of the throughput of each 
will be developed to compare the performance of each protocol. From this analysis, 
it is found that different fixed frame lengths provide better performance for various 
bit error rates. Although one frame length may be efficient at a given bit error rate, 
that same length may not be as efficient if the bit error rate changes. This is reflected 
in Figure 1.3 where, although a frame length of 1024 bits is most efficient at low 
bit error rates, at high bit error rates a substantial improvement in throughput is 


achieved using a frame length of 128 bits. 
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Figure 1.3: Dotted curves represent throughput of a fixed frame length ARQ proto- 
col for various frame lengths L. At low bit error rates, increased percentage overhead 
in small frames decreases throughput. At higher bit error rates. longer frames have a 
higher probability of error. The solid curve represents the desired performance from 
an adaptive ARQ scheme. 


To help improve performance through the entire range of bit error rates, an 


adaptive ARQ protocol can be implemented. With adaptive frame lengths. long 


we 





frames are used at lower bit error rates and shorter frames at higher bit error rates. 
By adapting to the appropriate frame length, the throughput achieved can approach 
the maximum at any BER of the channel (the solid curve in Figure 1.3). The concept 
of an adaptive ARQ strategy is proposed in [2]. This thesis presents the results from 
simulations where the length of a frame can adapt to the experienced bit error rate 


of the channel. 
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Il. RECEIVER ANALYSIS 


An analysis similar to that presented in [1] is carried out where, in an M-ary FSK 
communications scheme, one of M signals is sent during the time interval 0 < t < T. 
Customarily, M = 2* so that k-bits of digital information are represented by each 


signal. The possible signals to be transmitted are of the form 


Sm(t) = /2E, cos 2a fnt (2.1) 


where E, is the signal energy (E, = i s? (t)dt) and fm is one of the M frequencies 


used to distinguish among signals where 


=14 sien) 1). (2.2) 
Here, A; is an integer representing the spacing (relative to 1/T) between each of the 
M signals. This signal is modulated with an appropriate carrier frequency f, which 
is removed at the receiver. 

The receiver used in this analysis is depicted in Figure 2.1. After removing the 
carrier frequency and low-pass filtering (with cutoff frequency W > fmar). the signal 
is then sampled at a rate N/T (N is the size of the FFT) so that the outputs of the 
FFT correspond to a frequency spacing of 1/T. With this sampling rate, only the 
bins with a spacing Ay from Equation 2.2 will be used in the decision stage. 

The input to the receiver is the transmitted signal corrupted by zero-mean, 
Gaussian-distributed white noise n(t) with power spectral density (p.s.d.) of N,/2. 
The signal may also experience a frequency offset of f’ caused by doppler effects, 
receiver oscillator drift or other sources. In addition, the signal experiences slow 


Rician fading causing the amplitude a to be a random variable following a Rician 


7 











Figure 2.1: With a FFT receiver, the incoming signal is sampled, a window is 
applied, and the magnitude of the FFT output is used to determine which signal is 
received. 


distribution. Rician-distributed random variables in signal analysis are characterized 


by the ratio of the power of the direct signal component A* to the power of the 


diffused (or faded) component 207 





where Jo(x) is the modified Bessel function 


2n 


Io(z) s a an 


n=0 


At the receiver’s input, 


r(t) = /2E,acos(2x(fe+ fm + f’)t] + n(t) 


and after removing the carrier 


F(t) = V2E,acos [27( fm + f’)t}] + n(t). 


To simplify the signal representation, the following substitutions are made: 


q = 14+A;(m—1) 


kia 


p> 


8 





so that Equation 2.6 becomes 


F(t) = \/2E,a cos ae + 7(t) (2.9) 


where € now denotes the fractional (of 1/T) frequency offset and g represents the 
multiple of 1/T of the sent signal. With a sampling rate of N/T, the discrete form 


N 


As seen in the receiver of Figure 2.1, the discrete signal is windowed using a 
discrete windowing function w(n) and the FFT is taken. The outputs of the FFT 


can be expressed in real (denoted Xp(k)) and imaginary X;(k) parts 




















pel Qrk 
Xp(k) & - r(n)w(n) cos : uh 
n=0 N 
= ee tk hi 
= w(n) 2E,a cos MTT IM os zits + n(n)w(n) cos | (211) 
n=0 - 
N-1 ork 
X1(k) & r(n)w(r) sin — 
n=0 N 
— 2 Qn 
- {rein 2a cos “a+ li (n)w(n) sin | 212) 
n=0 


Keeping in mind that the channel experiences slow Rician fading so that a is 
assumed to be a constant over the period of interest 0 < ¢ < T, and, because the 
received signal is being corrupted by random noise following a Gaussian distribution. 
Xpr(k) and X;(k) are Gaussian distributed as well, since they are derived from linear 
combinations of Gaussian processes. As a result, the output of the FFT generator is 


a complex Gaussian process of the form 


X(k) = Xp(k) + pX7(k). (2.13) 





The necessary variable for the decision criteria will be the magnitude of the 
outputs of the FFT 

£ 1X (k)| = V/XR(b) + X7(4). (2.14) 

This is the envelope of the the complex signal X(k). As stated in [3], for a random 


complex signal X(k), the distribution of the envelope is 





d; 1 dk pt : 
fla) = Seesp |—sto (dt +12)) bo (SE) (2.18) 
where (overbar notation denotes expected value) 
w= [X(K)| = VV Xn(k) + X7(k) (2.16) 
and 
2 1 7 ——)? 7 
oh=5E (xt ~ X(b)| (2.17) 


The mean of X(k) (from Equations 2.11 and 2.12) for a fixed amplitude a are 
easily found by observing that n(n) is the only random quantity and is zero-mean. 


causing the second term of the sum to go to zero leaving the first term of the sum, 














N-1 
ba wed 2 ‘ k 
XR = ay2F, a oe Pra bad enn 
n=0 N N 
S mp (2.18) 
N-1 
oa : 2 cnk 
X; = ay/2k, » w(n) sin ay santa+ 6) sin n= 
n=0 
Bais (2.19) 
The variance in Equation 2.17 then becomes 
] N-1 eek 2 
oF = 5F | >> w(n)n(n) ee ae om n)sin ed | (2.20) 
n=0 n=0 








which, because of white noise, the autocorrelation function of n(n) is given by 
No 
R,(n, p) = Won p) (2.21) 
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where, again, W’ is the cut-off frequency of the low-pass filter at the receiver and N,/2 


is the p.s.d. of the additive white Gaussian noise. Equation 2.20 is reduced to 








WN s Ink Na _ o2ank 
a = = w*(n)cos? + w?(n)sin? —— 
a 2 py N dX N 
WN, N-~1 
ty >> w(n). (2.22) 
n=0 


Keep in mind that Equations 2.18 and 2.19 are still functions for a given k, m 
and ¢«. By factoring out a from mr and my, given in Equations 2.18 and 2.19, and 


substituting into Equation 2.16, 


m?, + mi 


a= 
II 


a’(mpg +m;) 


I> 


a’ By. (2.23) 


With slow Rician fading. Equation 2.15 is actually a conditional density on a and 





becomes 
d, 1 ad, Yk ‘ 
fa, {a( ds |@) — oper [sete + a) Io ( a : (2.24) 


The density function of d, can then be found by integrating over the Rician- 


distributed random variable a 





fa(dx) = f° fax(dsla) fa(a)da 
co dp -——(d24078 d ~ as (a? 4.47} : 
= i ae tot | a+ aA (: | ae 207 i (=) da 
0 ox ox oF OF 


a 1(@ a 
axa; ct ig ok oF 


ee ay h- + 1 ) ; 
x I ae (3 a to (SANE) 1 (24) a. (2.25) 


1 Bot + 0% 
6 ei ae ee CLL 8 ie (2.26) 
20% = 20% 20404 





and using the integral (from [4]) 


oo ] a? + b? ab 
~Q7 = ——exp | -———_. es 2.97 
[ re J,(az)J,(br)dz 5 7exP | 4Q2 | de (; i) ° (2:27) 


gives 





d, 1 dio}; + A®a% 
falas) = (<ieq) (ag) 0 |-heres 


x exp |( + “) : Ip SS) ; (2.28) 


ot of 4Q? 205,070)? 
From Equation 2.4, it is seen that Jo(z) is an even function. Replacing the Q? as 


defined in Equation 2.26 into Equation 2.28 results in 


Pdi, Sh nV ops ee, See ah 
dy Bro} + 0% 205,07 20%.03( 3,07 + 0%) 
ae area} (2.29) 
Brot + of 


which simplifies to 


(4 1 (2+ A? d, AVE ; 
fa, (dk) = (sorta) exp E (is ro Io Bo? + 0% : (2.30) 


Note that by making the substitutions 


Oo, = Boe tox 
AV Br, (2.31) 


Equation 2.30 is seen to be Rician-distributed since 


d,, 1d+az dpa 
fex(de) = exp (-3 a ) 1o( ig (2.32) 


k k k 


I> 


ak 








A similar result was found in [1] where the parameters o? and a, are now altered as 
a result of Rician fading. Note that if no fading is present (207 = O and A = 1), then 


Equation 2.32 is the solution found in [1]. 
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Ill. ERROR ANALYSIS 


Now having an expression for the distribution of the magnitudes of the FFT 
output, an error occurs when the desired frequency component (located at bin q of 
the FFT) is less than any of the other FFT bins of interest. With a sampling rate of 
N/T, the bins of interest are those corresponding to a frequency integer multiple of 
Ay. In addition, the FFT displays symmetry about the N/2 frequency bin so only 
the first N/2 outputs are of consequence. Finally, the first bin corresponding to a 
d.c. component is disregarded, as well as the bin located at N/2, since a trequency 
component in this bin will not satisfy the Nyquist criteria. 

As further discussed in [1] (referencing [5]), the union bound solution is carried 
out for the probability P, of symbol error. From Equations 2.18 and 2.19, the distri- 
bution of d, will differ with a positive or negative frequency offset ¢. It is assumed 
that a positive or negative frequency offset is equally likely. In addition. each of the 
M signals are assumed equally likely to occur with probability 1/4%. An error occurs 
when the magnitude d, of a FFT bin is greater than the magnitude d, (refering to 


Equation 2.7) of the bin of interest. The probability of symbol error is then given by 


1 
| ates OM De {Pr[d(q) < d(k)| + €,q] + Pr{d(q) < d(k)| — «.q}} 


k,q€1+ O¢(m — 1) 
m=1,2, ...M 
k #4 


where, from [5]. for Rician-distributed random variables 


(3.1) 


Pr(d(q) < a(k)} = 
Q( Va, Vb) — =——exp (-*) Io Vab) (3.2) 


c+) 








with the following definitions made 








a ay 
a= ; a 
a5 +o; 
2 
A ao 
Dos o? +0? 
qg k 
2 
A 9%, 
c= —2° 
o7. 


Evaluation of the Marcum-Q function is described in Apperdix A. Also, the depen- 
dency of €, g and k is found in the 6, term in the definitions made in Equation 2.23 
(with reference to Equations 2.18 and 2.19). 

Parameters for analysis include the direct-to-fading ratio 


a A 
DTF= 302 (3.3) 


and the signal-to-noise ratio at the receiver 


E A? + 20? 
CNRS a. 3. 
R N, N, (3.4) 


By solving each of the above equations for A?. it is found that 


o? SNR 4 


N, XDTF+1)? ee 
which is the diffused signal-to-noise density ratio. Using this, making the substitutions 


from Equation 2.31, and normalizing the transmitted energy (£, = 1 in Equation 2.18 


and 2.19), then a, 6 and c can be expressed in terms of the given parameters 


as 2DT FB, (3.6) 
= Be + Bg +77 W SN} w2(n) , 
2DTF 
b = Ba (3.7) 


By + By + W So w2(n) 


23,7 + WONT w?(n) 
2Bky + WSS wan) 
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A bound on bit error is then found by 


p< MP 





P,. (3.9) 


~ M-1 





IV. WINDOWING RESULTS 


Graphs of the probability of bit error as a function of the received signal-to- 
noise ratio (per bit) were constructed for various direct-to-fading ratios 4/20} for 
the rectangular and Hamming windowing cases. These results were generated from 
Equations 3.1 and 3.2 with evaluation of the Marcum-Q function and the modified 
Bessel function described in Appendix A. A plot with low fading channel (A?/20? = 
100) is given in Figure 4.1 below. The case studied in [1] (A; = 2. M = 8) is used so 
that a direct comparison of results can be made. In the case of low fading. the results 


are comparable. 
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Figure 4.1: Probability of bit error for low fading channel (a) rectangular. and (b) 
Hamming windowing. 32-point FFT, M = 8, A; =2 
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From Figure 4.1 the advantages of the non-uniform windowing can be seen as 


the frequency offset becomes larger for a fixed f,, and T. In both cases, as € grows 
linearly, the energy necessary to maintain a given bit error rate grows exponentially. 
However. the advantage of the Hamming window is reflected in a smaller exponential 
increase than the rectangular window. With lower values of frequency offset the 
signal energy necessary is larger than with uniform windowing, but, as € increases. 
the advantages of non-uniform windowing become apparent. 

A visual represenation of this was developed in [1] by graphing the amount of 
SNR degradation necessary to maintain a fixed probability of bit error as a function 
of increasing ¢ and is repeated here. Figure 4.2 gives the result for the low fading 


case with a fixed BER of 10-? (from Figure 4.1). 





5. = z ee Si 
/ i) 
' 
/ 
45+ F< 
/ 
| / 
a / 4 
i 
3.5; / 4 
_ : “A 
7 | 
= i. | 
3 3 Hamming ie | 
n . Le 
mm 25- ft 4 
‘al a. 
2 2 sere 3 / i 
& Le ! 
1.5 hee 7 = 4 
’ Ss 
; Rectangular i 








0.8: ‘i 
5 Direct,to-Fading Ratio = 100.0 | 
0 O01 02 03) 0405 


Fractional Frequency Offset (¢) 


Figure 4.2: Trade-off of uniform and non-uniform windowing for bit error rate of 
10-?. Non-uniform windowing becomes more advantageous as € increases. 
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Similar graphs of bit error rate as a function of SNR are given for the case 


of a slightly diffused channel (A?/20? = 3.16) in Figure 4.3. The results display- 
ing SNR degradation are given as well. From Figure 4.3 we see that with a more 
diffused signal, the cross-over fractional frequency offset decreases. This phenomena 
was studied closer by generating additional results in the same manner for various 
direct-to-diffused ratios and plotting the frequency offset cross-over point of each in 
Figure 4.4. From this graph we see that the jusfication to use non-uniform window- 
ing is very sensitive in the area of equal direct and diffused components. However, 
as the direct component approaches zero (Rayleigh fading), non-uniform windowing 


becomes more advantageous when minimal frequency offset is expected. 


Bit Error Rate 
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V. ARQ PROTOCOLS AND THEIR 
PERFORMANCE 


As seen in the previous section, noise on a communication channel can cause 
errors in the reception of signals. When an error occurs in computer communications, 
the receiving DLC must issue a command to the sender for retransmission of the frame. 
This negative acknowledgement (NAK) of the received frame, or acknowledgement 
(ACK) if no errors occured, can be a relatively short frame with a minimal chance 
of error or imbedded in a longer information frame with a higher probability of error 
occuring in the frame. As a safeguard, the original sending node assigns a maximum 
allowable waiting time, or timeout, where, upon its expiration, the frame is assumed 
to have been received in error and is then retransmitted. 

The efficiency of ARQ schemes is dependent upon the complexity of the DLC to 
implement it; the more complex the implementation, the higher the throughput rate. 
Current ARQ protocols—stop-and-wait, go-back-N and selective-repeat—are based 
on a fixed frame length. As will be seen from the throughput curves generated for 
each protocol, a higher performance may be achieved in areas of different bit error 


rates with different frame lengths. 


A. STOP-AND-WAIT ARQ 

The stop-and-wait ARQ protocol is the simplest protocol to implement and 
serves as a good example to illustrate ARQ schemes. In the stop-and-wait protocol, 
each frame is transmitted and the sending DLC then stands idle waiting for a response. 
If a negative response is made (either by a NAK or by a timeout). the frame is 


retransmitted. This process continues until all frames have been transmitted. As an 
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example, a communication is illustrated in Figure 5.1 between a sending DLC (node 
A) and a receiving DLC (node B). Here, frame ‘0’ is sent and received correctly by 
node B who ACKs the packet by requesting the next frame (frame ‘]’). An error was 
encounterd with frame ‘1’ and node B NAKs (negatively acknowiedges) the frame 
by re-requesting frame ‘1’ which is retransmitted by node A. The ACK from node B 
encounters an error and at the conclusion of the timeout, frame A retransmits the 
frame. 


{tprop fe fo-— timeout —_ot 


"delay 





ACK 
node B 


>| tack be 


Figure 5.1: Stop-and-wait protocol. Frame ‘0’ is sent and ACKed by node B. An 
error occurs with frame ‘1’ and is NAKed by B. The ACK from node B encounters 
an error, and at the conclusion of the timeout, frame A retransmits the frame. 


1. Performance Analysis Background 
Performance for each of the three protocols will be developed in their re- 
spective sections. Before proceeding with the analysis for the stop-and-wait protocol. 


the following background used in al] three performance analyses is developed. 


i) 
bo 








a. Assumptions 


The following assumptions are made throughout the performance analysis: 


e Node (A) continuously transmits to node (B). — With node A continuously 
transmitting, there are no queuing delays. Because node B has no information 


to send to node A, acknowledgement frames will only contain the ACK or NAK. 


e An acknowledgement is always received and is never received in error. — This 
assumption can be made because of the relatively short length of the acknow]- 


edgement frame. 


e The propagation delay is known. — This stipulation is only necessary to es- 
tablish an appropriate window size with the go-back-N and selective repeat 


protocols. 


e Processing time at each node is negligible and can be considered included in 


the propagation delay. 
e Bit error is independently and identically distributed for all bits within a frame. 


b. Average Transmission Attempts 
If the probability P, of bit error is independent and identically distributed for 
each bit within a frame (L bits long), the probability of a frame being received in 
error is given by 


Peal Pies 1 1 =P)”. (5.1) 


The probability that it will take 7 attempts to successfully transmit a packet can be 


expressed as 


Pr{number of attempts = 7] = P!"1(1— P.). (5.2) 





The expected number of attempts can be computed by ' 





E|[P(1-P.)| = (1- PE [Pi] (5.3) 
= (— PB) SP (5.4) 
t=1 
1 
a are (5.5) 


c. Throughput Definition 
Throughput is defined as 


a total bits successfully transmitted in a given time interval 


p (5.6) 


bit capacity of the channel in the same interval 
This can also be expressed as the ratio of the time to transmit a frame tian, to the 
average time to transmit a frame without error f. 

Because the information packet is appended with overhead bits £, some of which 
contains coding information, the throughput is scaled by this ‘coding’ rate. This 


results in a throughput efficiency given by 


4a ttrans (L = é) oie 
as Ll (5.7) 


2. Stop-and-Wait Performance 
The total time the transmitter expends for each frame is a sum of the trans- 
mission time of the frame, transmission time of the acknowledgement ft... and two 


propagation delays t,o) (which are assumed to contain the processing delay) for the 


round-trip cycle. This gives a total time for each transimission attempt of 


Te = daha 2b prag te tack: (5.8) 


The infinite sum is shown to be convergent by taking the derivative of the geometric infinite 
series solution where z < 1 


as, 91 ee 
bz 2? eae ~ (l-2z)? 











Using the results of Equation 5.5 for the average number of attempts, the average 


time to transmit a frame succussfully is given by 


FS trans + Lhorep + tack 


= 5.9) 
sw 1 Zs P. ( 9) 


By applying the definition of throughput given in Equation 5.7 and defining the 


number of frames that can be sent in one round-trip propagation delay as 


4 Lt ocp + tuck 


a : (5.10) 


ttrans 
the throughput of the stop-and-wait strategy is 


USSR EHO. Oa By ie!) 


Bee eas pl eg. 








(5.11) 


Theoretical curves for Equation 5.11 are given in Figure 5.2 for various frame lengths. 
The simulations described in Appendix B were performed and the results are given 
here as well. 

Although the stop-and-wait protocol is easy to implement, the efficiency of the 
system is compromised by waiting for a response after each packet is transmitted. The 
following two protocols do not stand idle waiting for a response and, subsequently, 


have better expected throughput than the stop-and-wait. 


B. GO-BACK-N 

In the go-back-N protocol, node A continuously sends frames to node B. Node 
B keeps node A aware of its status by responding with an acknowledgement which 
contains an RN (received number) indicating that all frames preceding RN have been 
received correctly. RN also indicates that all frames from RN onward need to be sent 
by A. As frames are received without error, RN is incremented. 

In the event of an error, RN will remain constant as illustrated in Figure 5.3. 


There is a limit to the number of frames node A can send beyond RN so that a frame 
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Figure 5.2: Performance of stop-and-wait for frame lengths of 1024. 512. 256 and 
128. Channel capacity is 32000 b/s and a round-trip propagation delay of 0.02 sec- 
onds. 


in error can be detected and retransmitted by node A. The number of frames that 
can be transmitted beyond RN is the Nin go-back- N and is referred to as the window 
length. Because RN remains constant, when an error occurs, node A will stop as soon 
as the end of the window is reached (SV +1 = RN +N) and go back N frames 
and begin retransmission of frame RN and all subsequent frames. To minimize the 
number of frames needed to be retransmitted in the event of an error, the value of .V 
is set to the number of frames that can be transmitted in the time interval given in 


Equation 5.8. Using Equation 5.10, 
Nz=a+tl. (5.12) 


With the go-back-.V protocol, the transmitter can send frames continuously. 
However, upon each frame received in error, the protocol requires that the sending 


DLC go back .V frames and resend them. Thus, when successful on the ath attempt, 


26 





Go-back-4 


oli j2tsiats[2tatalsio[z| 


_ 


node B 






Figure 5.3: Go-back-N protocol for N = 4. With frame 2 received in error. RN 
stays constant at 2 until node A has reached the end of the window. Node A detects 
the error and must "go back 4” frames. 


the DLC has sent a total of 1+ (7 —1)N frames. The average number of attempts 
was given in Equation 5.5 which results in an average number of frames (f) that need 
to be sent of 


1 
Pest evans 
a 1) 


1+ P(N —1) 
1-P, - 


SHI 
MI 


Using Equation 5.12, the average time to transmit a frame successfully is 


t bivanat 
Livanall ag aP.) - 
for (5.14) 


Applying the definition given by Equation 5.7, the theoretical throughput for the 


go-back- N protocol becomes 














where the dependency of P, on L is shown in Equation 5.1. A graph of theoretical 
throughputs for various frame lengths is given in Figure 5.4 with the simulation 
results superimposed. It can be seen from the graph that an adaptive strategy will 
not improve the performance of the go-back- N protocol. A frame length of 1024 could 
be used at all bit error rates and only a fractional percent of throughput would be 


lost at higher bit error rates. 
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Figure 5.4: Performance of go-back-N for various frame lengths. A frame length 
of 1024 could be used at all bit error rates with only minimal loss of performance at 


higher BER. 


C. SELECTIVE-REPEAT ARQ 

The selective-repeat protocol is similar to the go-back-N except that when an 
error occurs, only the frame that is in error needs to be retransmitted. Here it 
is important to acknowledge each frame (both positive and negative) so that the 
transmitter accounts for all frames sent out. If a response is not received. that frame 


reaches its timeout and is retransmitted. 
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node B 
Figure 5.5: Selective repeat protocol. 






The obvious advantage of selective-repeat over go-back-N is that only the frame 
in error is retransmitted as opposed to sending all following NM — 1 frames in the 
go-back-N protocel. However, this requires that the receiving DLC have adequate 
memory to store frames that it receives and to be able to sort them as they become 
available. If memory in the receiver is full, then all incoming frames are NAKed, 
regardless of any errors. 

The throughput analysis for selective-repeat (with an infinite receiver buffer 
size) is straightforward from Equation 5.5. Because the only delay experienced is 


that of the retransmission of a single frame when an error occurs, 


tirans 


t= 
1— P, 


(5.16) 


which results in a theoretical throughput for selective-repeat with an infinite receiver 
buffer as 


ff 
Psr = (1 Sa (5.17) 


Plots of theoretical throughputs with infinite receiver buffer length are given in Fig- 


ure 5.6 with the results obtained in the simulation. 
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Figure 5.6: Performance of selective repeat with an infinite buffer. 


By limiting the size of the buffer to the number of frames that can be transmitted 
in a round-trip delay N, a lower bound on the throughput is given in [6] by 


Xo 


> ——_——__ 5.18 
PEN ae hee fad (5 ) 
where 
B N-1 
MG oe Ba) 
1 = Pek (ake? xe Ppt? 4: Pix?) 
do = 3-— P? aN? _ a? 3? = pry? 
and 
a= 1— P? 
Bo LaP? 
> ee 
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Figure 5.7: Performance of SR with finite receiver buffer (fixed length L = 1024). 
N = 22 is the number of frames transmitted in a round-trip delay. From the analysis 
of [6], a lower bound on throughput is given. 

A comparison of the upper bound (infinite receiver buffer) and the lower bound given 
in Equation 5.18 is given in Figure 5.7 for a fixed length of L = 1024. For the 
simulation, a receiver buffer size of 64K is used instead of restricting it to .V (.V = 22 


for the parameters of the simulation). A buffer size of 64K was believed not to be 


too excessive, yet imposes some restriction on the system. 
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VI. ADAPTIVE ARQ PROTOCOLS 


Implementation of the previously described ARQ schemes remains the same in 
an adaptive environment. The adaptation takes the form of varying the length L 
of the frames. With an increased bit error rate, the increase in frame error rate in 
Equation 5.1 can be offset by decreasing the length of the frame. Alternatively, if 
the bit error rate of the channel decreases, a longer frame length will decrease the 
percentage of a frame dedicated to overhead (assuming the overhead remains constant 


for varying lengths). 


A. ADAPTIVE STOP-AND-WAIT USING A SEMAPHORE 

The decision to change the length of the frame was made in three different 
schemes. The first scheme is taken from [7] and invokes a counting semaphore m 
that is incremented each time a frame is acknowledged and decreased each time a 


frame needs to be retransmitted. The length of the frame is then dependent upon 


the current value of the semaphore by 
Ll yw>wsm<n 
204% Sm< 42 
Eee (6.1) 
nl Yn-1 < ™m = Yn- 

For this thesis, /= 128 and n = 4. 

Figure 6.1 shows the gain associated with an adaptive stop-and-wait protocol. 
This scheme adds little complexity to the DLC implementation and provides the 


desired results. 
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Figure 6.1: Adaptive stop-and-wait performance. Improved performance is obtained 


in the region of high BER. 


B. ADAPTIVE SELECTIVE-REPEAT USING THRESHOLDS 


Another criteria compared the number r of frames transmitted after a fixed 
value 6 of frames were received in error to a threshold to determine if the frame length 
should be increased (doubled) or decreased (halfed). Two schemes were implemented 


where either a single threshold 


L/2 b/r>y7 
eo { 2L bfr<y7 (6.2) 


or a dual threshold 
L/2 prj > pr 
L= 4 2h pa > pr (6.3) 
L otherwise 


was used. The computation of pz 2. pr and p21 are based on the experienced bit error 
rate. The experienced P, can be computed from the ratio of transmitted frames to 
frames received in error by applying the inverse of Equation 5.1] 


on ae 





1-P,=(1-P,)' 


ere = 
r 
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1/L 
A=1-(=] (6.4) 


Throughputs for lengths 2L and L/2 can then be computed for the selective repeat 


protocol with Equation 5.17. The decision criteria of Equation 6.3 then simplifies to 


L/2 1-€/(L—2) > ((r —b) fr}? 
ff | QL 1-e/(2L—2) <(r—b)/r (6.5) 
L otherwise 


Table 6.1 can be constructed for the dual threshold selective-repeat case. The deter- 
mination of y in Equation 6.2 was taken as an average (7 = 0.88) of the thresholds 


necessary to cause transition from table 6.1. 





TABLE 6.1: Adaptive SR transition levels. The necessary value of (r — b)/r to 
cause transition either up or down with the dual threshold implementation (¢ = 32 
bits). 


The performance of the adaptive selective-repeat for a fixed buffer length is 
given in Figure 6.2. Implementation of this scheme further increases the complexity 
of the selective-repeat protocol. Now the transmitting node must keep track of the 
length of each frame it is sending out. In addition, re-arrangement of the frames at 
the receiving node becomes quite complex. For example, if a frame is NAKed and 
the decision to change frame lengths is made, the packet sent in error will now be 
divided between two frames each of which may or may not be received correctly in 


order upon retransmission. 
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Figure 6.2: Adaptive selective-repeat. Both dual and single threshold results are 
given. Performance of the single threshold implementation is slightly lower in the 
areas where the decision to change lengths can “toggle” between two lengths. 


C. GO-BACK-N ARQ 

The go-back-N protocol does not lend itself to an adaptive environment as was 
seen in Figure 5.3. The difference in throughput rates between the largest frame 
length used (1024) and the maximum achievable throughput is negligible. This is 
because the inefficiency of go-back-N arises from the necessity to go back N’ frames 
when an error in a frame occurs. This is depicted in Figure 6.3 where, if the first 
frame is in error, the N — 1 following frames composed of L bits each result in a 
possibly unnecessary retransmission of (N — 1)L bits. If the frame length is halved 
(and N doubled to satisfy Equation 5.12), the probability of having to re-transmit 
the same number of bits is the same because the probability of transmitting two 
successive frames of length L/2 correctly is equal to the probability of transmitting 


a single frame of length L without error. Although some advantage may be gained 





in successfully transmitting the first frame of length L/2, there is also an increase 
in percentage overhead for each frame. For example, with a BER of 1074, an initial 
length of L = 1024 and overhead ¢ = 32, if the frame length is halved, Psucc increases 
by only 4.7% while the percentage of frame overhead increases 3.1%. This results in 


negligible improvement in the areas of higher bit error rates. 


Repeated if preceeding frame in error 





Figure 6.3: Inefficiency of the go-back-N protocol occurs because frames following an 
error must be repeated. If the frame length is halved, the probability of retransmitting 
the same number of bits remains the same. 


36 








VII. CONCLUSIONS 


In an M-ary FSK scheme implemented using sampled data at the receiver and a 
FFT to determine which signal was sent, it was found that non-uniform windowing can 
be used to increase the performance of the system in the presence of large frequency 
offsets. Implementing non-uniform windowing causes greater separation between fre- 
quency components of the FFT when frequency shifts occur between transmission 
and reception. Non-uniform windowing also causes attentuation of the signal making 
detectability more difficult. A balance between improved frequency separation and 
attenuation is achieved as the frequency offset becomes more dramatic. In the case 
of a fading channel, the frequency offset necessary to justify non-uniform windowing 
becomes smaller as the channel approaches Rayleigh fading. 

The application of an adaptive ARQ protocol was shown to improve the through- 
put of the stop-and-wait and selective-repeat ARQ protocols. With an adaptive frame 
length, throughput can be increased by decreasing the probability of frame error in 
high bit error rate situations through a reduction in the frame length. With low 
bit error rates, a larger frame length reduces the percentage overhead in the frame 
thereby increasing the throughput. The go-back-N protocol does not lend itself to 
adaption. From theoretical analysis and simulation, it was found that a fixed frame 
length can be used in the go-back-N protocol through all bit error rates with minimal 


loss of performance. 
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APPENDIX A: ERROR COMPUTATION 


The probability given in Equation 3.2 was computed using the enclosed C pro- 
gram code. Values for the Marcum-Q function and the Modified Bessel Function were 
calculated to within a designated tolerance (10-°) using power series expansions of 


the functions from [4] and [5] 


] if 3=0 
eur ifa=0 

Q(a, 8B) = (1 + Ip(a?))e-™” ifa=23 (A.1) 
exp |—2+2 








fo} 
n= 
fore) 
n=l 


lL, 1 if2= 0 Aa 
E)= (¢/2)" 52 eer otherwise (A.2) 


"In(ad) ifa<B 
"In(aB) ifa>3Z 


where 
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1. Bit Error Program Code 


[PRREEEREERAREERRESERERESESEREEESEREEEEERER GEER EEREREEE REESE 
* Pefad.c 

Calculates the probability of bit error for a Rician 

fading channel and a DFT based receiver for MFSK 


++ &* & He 


by: Chris G. Kmiecik 

date: 22 April, 1990 * 
EREEREREEEEERERREREREREREREREREREREEEES EREEEEEEERER SE EERE ES / 

#include <stdio.h> 

#include <math.h> 


+ + & & 


#define pi 3.141592654 
#define sqrt2 1.414213562 
#define epsilon 0.00000001 


#define N 32 
#define M 8 
#define 1 3 
#define W 16 
#define df 2 


double window([N]; 


main() 

{ 

double mip(), mrp(), I(), Q(), power(), rectangle(), hamming(), hanning(); 
double Perror, lastPerror, Enoise; 

double a, b, c, numb, numc, denon; 

double meanrq, meaniq, betaq; 

double meanrk, meanik, betak; 


long fact(); 
FILE *output; 
int n, KX, q, einc; 


float SNRdb, SNR, DTF, DTFinc, ratio, e, freqe; 


system("rm PefadoutHAMM.mat"); 


Enoise=0.0; /* «/ 
for (n=O ; n <= N-1 ; n++) { /* Compute the values */ 
window(n+1]=hamming(n) ; /* for the window at */ 


Enoise=Enoise+Wewindow[n+1]*window[n+1]/2.0; /* each ’n’ and the energy */ 


} 








for (DTFinc=2.0 ; DTFine >= -2.0; DTFine = DTFinc-0. 


DTF=pow(10.0,DTFinc) ; 
for (e=0.0 ; e <= 0.5 ; e=et0.1) f 


Perror=3.0; lastPerror=4.0; 
SERdb=2; 


while((Perror >= 0.00001) 
&&(SERdb <= 99) 
&&(lastPerror-Perror >= .00001)){ 
lastPerror=Perror ; 
SNRdb++; 
SNR=1*pow(10.0,SNRdb/10.0); 


ratio=SNR/(DTF+1) ; 
Perror=0.0; 


for (q=1 ; q <= M*df-1 ; q=qtdf) f{ 
meanrg=mrp(q,e,q); 
meaniq=mip(q,e,q); 
betaq=meanrq*meanrq+meaniq¢*meanigq; 
numc=0.5*betaq*ratio+Enoise; 
numb=(SNR-ratio)*betaq; 


for (eine = 1; eine <= 2; einc++) { 
freqe=etpower(~-1.0,einc); 


for (k=1 ; k <= (N/2)-1 ; k=k+dt) { 
if (k !=q) { 


/* 


/* 


/* 
/* 
/* 
/* 


/* 


/* 
/* 
/* 
/* 
/* 
/* 


/* 
/* 


/* fading ratio 


Cycle through various 
frequency shifts 


Continue computing Pb 
until SNR becomes large 
or error becomes small 
or error does not change 


Convert Eb to Es 


Each value of q donates 
to Pe. First 


25){/* Specify Direct to */ 


*/ 


*/ 
*/ 


*/ 
«/ 
*/ 
*/ 


*/ 


*/ 
*/ 


find the mean of the real*/ 
and imaginary parts then */ 


specify all parameters 
changing with q only. 


Include both plus and 
minus Doppler shift 


meanrk=mrp(q,freqe,x); /* Calculate values in 
meanik=mip(q,frege,k); /* equations 2.17, 2.18, 
betak=meanrk*meanrkt+tmeanik*meanik; /* 2.22, 3.1. 


c=numc/(0.5*betak*ratio+Enoise) ; 
denom=0 .5*(betaqtbetak) *ratio+2.0*Enoise; 


a=(SNR-ratio)*betak/denom; 
b=numb/denom; 


Perror=Perror+Q(sqrt (a) ,sqrt(b)) 
-(c*exp(-0.5*(atb))#*1(0,sqrt(a*b)))/(c+1); 


} /* if (k != q) */ 


/* Equation 3.2 


} /* for (k=1 ; k <= (N/2)-1 ; k=k+df) */ 
} /* for (eine = 1; einc <= 2; einct+) */ 
} /* for (q=1 ; q <= M*df-1 ; q=qtdf) */ 


Perror=Perror/(2.0#M); 
if (Perror <= 0.2) { 


output=fopen("PefadoutHAMM.mat","a") ; 
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/* Equation 3.1 


*/ 
*/ 


*/ 
*/ 


*/ 
*/ 
*/ 


«/ 


*/ 








fprintf(output,"%10.8f %4.1f %4.2f %5.2f\n",Perror,SNRdb,e,DTF) ; 
fclose(output); 
} 


} /* while (Perror >= 0.00001) */ 
} /* for (e=0.0; e <= 0.5 ; e=e+0.1) #/ 
} /* for (DTFinc=2 ; DTFine >= -2; DTFine = DTFinc + 0.02) */ 
} /* close main */ 


4] 














[ERERREEAEEERAREEEEERERERSERESEERE EE EERE SERRE EERE SEE EEE 


* * 
* Function mrp * 
* * 


ERERERAREREEERESEESEEA SEES ESE ESEEEEEEEESEEESEEE EEE EE / 
double mrp(q,e,k) 


int k,q; 
float e; 

{ 

int n; 


double ans; 

ans=0.0; 

for (n=0; n <= N-1; n++) { 
ans=ans+cos(2*pi*n*(qtetk)/N)*window([n+1]+cos(2*pi*n+(qte-k)/N)*window([n+1] ; 
} 

if ((e == 0.0) && (k != q)) ans = 0.0; 


ans=ans/sqrt2; 


return(ans); 

[EERE REERAEREAREEEEEEREER EKER EER KEE EAREKEE EEE REE EE 
* * 
* Function mip * 
* * 


EKER EER DEER RE EE ER a a / 
double mip(q,e,k) 


int k,q; 
float e; 

{ 

int n; 


double ans; 


ans=0.0; 
for (n=0; n <= N-1; n++) { 
ans=ans+sin(2*pi*n*(qt+te+k) /N)*window(n+1]~sin(2*pi+*n*(qte-k)/N)*window([n+1]; 


} 
if ((e == 0.0) && (k != q)) ans = 0.0; 
ans=ans/sqrt2; 


return(ans); 
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[RERREREEREEREEEEEEEEREREEEEEEEEECE SESE EEE REE EEESESES 


* * 
* Marcum-Q Function * 
* * 
ERERE EERE EERE REESE REE EER EERE EE EER EEE EE EE ER ES EEE EEE / 

double Q(alpha, beta) 

double alpha, beta; 

{ 

double ans, e, inc; 

int n; 


if (beta == 0.0) 
ans = 1.0; 
else if (alpha == 0.0) 
ans = exp(-0.5*beta*beta) ; 
else if (alpha == beta) 
ans = 0.5*(1.0+I(0,alpha*alpha)*exp(-1+*alpha*alpha) ) ; 


else if (alpha <= beta) { 
ans = 0.0; 
e = exp(-0.5*((alpha*alpha)+(beta*beta))); 
n=0; inc = 1.0; 
while (inc >= epsilon) { 
inc=e*power (alpha/beta,n)*I(n,alpha*beta) ; 
ans = ans + inc; 


ntt; 
} 
} 
else { 
ans = 0.0; 


e = exp(-0.5*((alpha*alpha)+(beta*beta))); 
n=1; inc = 1.0; 
while (inc >= epsilon) { 
inc=e*power (alpha/beta,n)*I(n,alpha*beta) ; 
ans = ans + inc; 
n++; 
} 
ans=1.0-ans; 


} 


return(ans); 


} 
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[RRRERAEREREESEEAEEREREREEESEAAEREREREREREEERREEEEEES 


* 
* 
* 


double 
int 
double 


* 
Modified Bessel Function * 
* 
CERERREEKEEREEREREREREEEEEREE EERE EERE SEEEREEEREE EERE / 
I(v,x) 
Vv; 
x; 
ans, inc; 
n; 


else { 
x=x/2 


.0; 


ans=power(x,v)/fact((long) v); 


inc=ans; n=1; 


while 
inc 
ans 


F 
} 


(ine >= epsilon) { 


=incex*x/(n*(v+n)); 
=anstinc; 
nt+t; 


» 


return(ans); 


} 
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[ERERREEREREREEEEREEEERAEEERAEAEEREEEE REESE EEEERERERE EE 


* * 
* Function rectangular * 
* * 


Etiiit iter i iii iii tere oe ri er 
double rectangle(n) 


int n; 

return(1.0); 

} 

[RERERREREREREEREERERESEEEREEE EERE EERE EEREE EEE EEEERERE 
* * 
* Function Hamming * 
* * 


WAR RE RR RE OE EER ER EE EE RE EE / 
double hamming (n) 
int n; 


{ 
return(0.54-0.46*cos(2*pi*n/(N-1))); 


e 








[RERERERERREEEREREEREEEEEEREEEEEEEREER EE EEE RERAEE EERE 


* * 
* Factorial * 
* * 


EEREERAAKEREREEEEREE EERE RERERRESEEEEEEEEEERERE EEE EEE / 
long fact(x) 


long x; 

{ 

long i, ans; 

ans=1; 

if (x == 0) 
ans = 1; 

else 


for (i = 1; i <= x; it+) ans = ans*i; 


return(ans); 

VALI TIES TSP PETS e TPP PPP PS SPP T PT PPPS PEP PPC Se SSP SSS Try 
* * 
* POWER * 
* * 


SILL I TT TT TTT Tir tT Tole TT Tt TT TTT TTT TTT TTT TT TT eT TT Tea 
double power(x,n) 


int n; 

double x: 

{ 

double p; 

int recip; 

if (n <0) { 
recip=1; 
n=abs(n); 
} 


for (p=1.0; n > 0; --n) P=p*x; 
if (recip == 1) p=1.0/p; 


return(p); 


} 
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APPENDIX B: ARQ SIMULATIONS 


The ARQ simulations were carried out using a C-program. Although one node 
would be transmitting to another in the application of an ARQ protocol, the simu- 
lation only involves one node. A doubly-linked list maintains a record of jobs to be 
performed in chronological order by the node. The nede can either transmit (SEND’) 
or receive an acknowledgement (’RECV’) for a given frame. 

For each transmission, the time of completion of transmission and the computed 
time of reception (based on Equation 5.8) are added to the list of jobs. At tne 
completion of transmission, the node then transmits the next frame if the protocol 
allows it. 

At the execution of "RECV’, a uniformily distributed, pseudorandom number 
between zero and one is compared to the theoretical frame error rate (Equation 5.1). 
If the comparison dictates, a NAK is recorded by the node and the appropriate 
retransmission protocol is applied. Particulars for implementing retransmission in 
each protocol are given below. 


1. Selective-Repeat 


This simulation served as the template for the others (see Section B.4) in that go- 
back-.V and stop-and-wait are “subsets” of this program . The adaptive nature of the 
program is contained in a section added after receiving each frame's acknowledgement. 


(see Figure B.1). 





2. Go-Back-N 


Before sending a frame, the window is checked to see if it is full. If the window 


is full SIN is reset to the beginning of the window (Figure B.2). 

3. Stop-and-Wait 

Here, only the time to receive each packet is loaded onto the linked list. When 
a frame is received, it is either ACKed or NAKed and a frame is sent out (see Fig- 
ure B.3). If the frame is ACKed, the number of successful bits is incremented. 


4. Selective-Repeat Program Code 


[RRRERERREEAAREKREEEKEEREKEERESEREREKEEEEKERERE EERE ERA RETA EE EEE EEE E EERE EE AES 
* sr.c 


Simulation of selective repeat protocol with data generated 
for BER from 1E-6 to 1E-2. 


Receiver buffer size is kept track of by the sender 


* &£ &£ & & & & 


by: Chris G. Kmiecik, LT, USCG 

date: 0S February, 1990 * 
SEI OOOO OOOO RR ao / 

#include <stdio.h> 

#include <malloc.h> 

#include <math.h> 


+ *£ eH HH He 


#define MALLOC(x) ((x *)malloc(sizeof(x))) 
#define max_FN 4096 
#define max_buffer 65536.0 


typedef struct job { 
float time, length; 
char type[5] ; 
int FN; 
struct job *previous, *following; 
} job_type; 


float LogBER ; 
job_type *top_job; 
FILE *xof, *tf, *rf, *ttf; 


main() 


{ 


framestatus [max_FN] ,framelength[max_FW] ; 
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int queue; 


int buffer, begin. window, SN, i; 

float completed, m, Ps; 

float BER, simulation_length; 

float capacity, propogation, overhead, ack_length; 
float rand, length; 

float successful_bits, throughput, tr, tt; 
job_type *temp, *previous_job, *job, *newjob(); 
unsigned long int seed; 

char *strcat(); 

seed=1234; 


system("rm throughsr"); 


simulation_length=250.0; /* 
capacity = 32000; /* Simulation Parameters 
propogation = 0.7; /* 
overhead = ack_length = 32.0; /* 


for (length =1024.0; length >= 128.0; length=length/2.0) f{ 
for (logBER = -6.0; logBER <= -2.0; logBER=logBER+.1) { 


BER = pow(10.0,1lo0gBER) ; /* 
begin_window = queue = SN = i =1; /* Initialize variables 
successful_bits = buffer = 0.0; /* 


while (i <= max_FN) framestatus[i++]=0; 
m=completed=0; 


top_job = newjob(0.0,"send",0,0.0); 
top_job~>previous = top_job; 
previous_job = top_job; 


job = newjob(simulation_length+1,"send",0,0.0); 


job->previous = previous_job; 
previous_job->following = job; 
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*/ 
*/ 
*/ 
*/ 


*/ 
*/ 
*/ 


[RERREEEEREERERERESEESE EES Begin Simulating EREREEEREEERERE ER EREEEER ES EE EE EE / 
while (top_job->time <= simulation_length){ 


[RRRKEREEEREREREREREEEEEEE EES ER EEE 
* * 
* Time to Transmit a Packet * 
* * 
REREEEEREREREERERERS ES EERERES SESE / 
if (!(stremp(top_job->type,"send"))){ 


i = begin_window; 
queue = SN; 
while (i != SN) f 
if (framestatus[i] == 3) { 
queue = i; 
break; 
} 
if (++i > max_FN) i=1; 


} 


/* * 
* Re-transmit Packet in Queue * 
* * 
EKRREREAKRERAAEEEREEEREERERERE EEE / 
if (queue != SN) { 
tt = length/capacity; 
job = newjob(top_job->time+tt ,"send",0,0.0); 
place_job(job); 
tr = propogation+((lengthtack_length)/capacity) ; 
job = newjob(top_job->time+tr,"recv", queue, length); 
place_job(job); 
framestatus [queue] = 4; 
framelength[queue] = length; 
} 
/* * 
* Transmit if Window not Full * 
* * 
Ae eR a a OE OE OE / 
else if (buffertlength < max_buffer) { 
tt = length/capacity; 
job = newjob(top_job->time+tt,"send",0,0.0); 
place_job(job) ; 
tr = propogationt((length+ack_length)/capacity) ; 
job = newjob(top_job->time+tr, "recv",SN, length) ; 
place_job(job) ; 
buffer=buffer+length; 
framestatus(SN]) = 1; 
framelength[SN] = length; 
if (++SN > max_FN) SW = 1; 
} 











/* * 
* IF Window is full * 
* Try to Send After Receiving * 
* Next Packet * 
* * 
FEREEKREEREEERERE SEES ESEE ES EERE EE / 
° else { 


job = newjob(top_job->following->time,"send",0,0.0); 
place_job(job); 
, } 


} /* Packet transmission complete */ 


[EREEREEEEREREEREERESEEREREREKEREE 


* * 
* Time to Recieve a Packet * 
* * 


Seto a Rk Ok a  / 
else if (!(stremp(top_job->type,"recv"))) { 
seed=seed*1103515245+12345; /* Psuedo-random +*/ 
rand=(seed/65536) % 32768; /* generator */ 


/* * 

* NAK the received packet * 

* * 
. SRA TERE EERE ER REESE EE / 


if (rand/32767 > pow((1.0-BER) ,top_job->length)) { 
framestatus[top_job->FN]=3; 


} 

/* * 
* ACK the received packet * 
* * 
SAE OE I a a ik / 
else { 

framestatus [top_job->FN]=2; 
} 


} /* Packet Received */ 
[REAR A HEE EE 
* * 
* Update Window Beginning * 
* * 
HEREREREERRERER ES EERE EEE EE ER ERE / 
while (framestatus[begin_window] == 2) { 


successful_bits=successful_bits+framelength [begin_window]~-overhead; 
framestatus [begin_window]=0; 
buffer=buffer-framelength[begin_window] ; 

if (++begin_window > max_FN) begin_window=1; 


} 





[/REKRRERKESEREEERERERERESSESEEES EEE 


* * 
* Get next job * 
* Free current job memory * 
* * 


ERAREREEEEREEREREREE EERE EEE EE ESE / 

temp=top_job->following; 

free(top_job) ; 

top_job=temp; 

} /* repeat for next job -- “while(top_job->time<=simulation_length)” */ 


/eeeeeEEEEERERRES ORES ESRESES Simulation Completed REET KEE EEE EE EEE EEE EER / 


[REREREREREREREASERERER EERE RERERE 


* * 
* Record Throughput Results * 
* * 


SE Eo oo et / 

ttf = fopen("throughsr","a"); 
throughput=successful_bits/(capacity*simulation_length) ; 
fprintf(ttf,"46.5f %11.8f\n", throughput, BER) ; 
fclose(ttf); 


} /* next BER -- "for (logBER=-4.0; logBER<-2.0; logBER=logBER+.1)}" */ 
} /* for (length =1024.0; length >= 128.0; length=length/2.0) */ 


} /* close main */ 








[RRREEREERREEEERAESE EER EE EEE EEEEE EERE EEEEE ESSERE EERE 


* * 
* Function newjob * 
* * 


HERRERA EKR EERE REESE EEEREEEEEEEEREEREEE SEE SEEKS / 
job_type *newjob(time, type, FH, length) 


int FN; 

float time, length; 
char type(5]; 

{ 


job_type *new; 


if (!(new=MALLOC(job_type))) { 
fprintf(ttf,"out of the storage \n"); 
fclose(ttf); 
fclose(rf); 
fclose(of); 
fclose(tf); 
exit(1); 
} 


new->time=time; 
strcpy(new->type,type) ; 
new->FN=FN; 

new~->length=length; 
new->previous=new->following=NULL; 
return (new) ; 


} 
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/SEKEKKEREAKE REE ESERSESEEAA REESE ESSERE ESE EEELE REESE 


* * 
* Function place_job * 
* * 


ERERKERSERESEEESESSESEREEERERESREREREREREEEEEEESEEEEE / 
place_job(job) 

job.type *job; 

{ 


job_type *current; 


current=top_job; 
while (job->time >= current->time) f{ 
current=current->following; 
if ((current->following == NULL) && (job->time >= current->time)) { 
current->following=job; 
job->previous=current ; 
job->following= NULL; 
return; 
} 
¥ 


current->previous->following = job; 
job->previous = current~>previous; 
current~->previous=job; 
job->following=current ; 


} 
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/* * 
* Decision Criteria for Level * 
* Change * 
* * 
SOO IGOR IO a oI ik aK i kak kak / 
if (m >= 1) { 
Ps=(completed-m) /completed; 
if (Ps < 0.87) { 
length=length/2.0; 
: Ciength < 128.0) length=128.0; 


else { 
length = length*2.0; 


: (length > 1024.0) length=1024.0; 


m=completed=0 ; 


} 


Figure B.1: Decision segment for single threshold. 








DR IIR IO IOI I IRR IO Ik 


* x 
* Time to Transmit a Packet * 
* * 


FE SASSO OSORIO IO OIG IE io / 
if ('(strcmp(top_job->type,"send"))){ 


/* * 
* Transmit if Window not Full * 
* * 


TTTrrttt ttre ttt tet Teter ret. Ste to) 
if (SN-RN <= window) { 


/* * 
* IF Window is full * 
* Go-back-N * 
* * 
SOOO IO OR OR IO OR II IORI ia ac / 
else { 

SN=RN=1; 


job = newjob(top_job->time, "send" ,0,0.0); 
place_job(job); 
} 


} /* Packet transmission complete */ 


Figure B.2: Go-back-N implementation. 
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ROI GIO OR IG I I IR aK 


* * 
* Time to Recieve a Packet * 
* * 


rVTTTTTTTTTttrTT tt rT Tit tT Tite tte ttt a 

if (!(stremp(top_job->type,"recv"))) { 
seed=seed*1103515245+12345; 
rand=(seed/65536) % 32768; 


/* * 
* ACK the received Packet * 
SAO IOI I IGRI OK IG i kak kak tok ak ak ak akc // 
if (rand/32767 < pow((1.0-BER) ,top_job->length)) { 
successful_bits=successful_bitsttop_job->length-overhead; 
if (SN == 1) SN=0; 
else SN=1; 


tr=((length+tack_length) /capacity)+propogation; 
job = newjob(top_job->time+tr,"recv",SN, length) ; 
place_job(job); 

} 


temp=top_job->following; 
free(top_job) ; 
top_job=temp; 

} /* repeat for next job */ 


Figure B.3: Stop-and-wait simulation block does not require a "Time to Transmit 
a Packet” section. 
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